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ABSTRACT 

The origin and the profile of the iron Ka Hne from active galactic nuclei is being stud- 
ied. The model with a quasi-spherical distribution of clouds around a compact galactic 
nucleus is explored in which the intrinsic spectral feature around 6-7 keV is substan- 
tially affected by Comptonization of primary X-rays. The predicted spectral profiles 
are influenced mainly by motion of the clouds in strong gravitational field of the centre, 
and by orientation of the source with respect to the observer. The case of the Seyfert 
galaxy MCG-6-30-15 is examined in more detail and free parameters of the model are 
determined. A robust result is obtained: acceptable fits to the observed profile corre- 
spond to the spherical distribution of the clouds with strong self-obscuration effects 
and with rather high value of ionization parameter. This is a different situation from 
that considered in the model of a cold illuminated disc. Nevertheless, geometrically 
flat (disc-type) configurations are not rejected. 

Key words: Galaxies: active - Galaxies: Seyfert - X-rays: galaxies - X-rays: stars - 
Accretion, accretion discs - Galaxies: individual (MCG-6-30-15) 



1 INTRODUCTION 

An ever-increasing accuracy of the determination of the ob- 
served shape of the iron Ka line in the Seyfert 1 galaxy 
MCG-6-30-15, and the evidence of the broad and skewed 
profile (Tanaka et al. 1995; Iwasawa et al. 1996) have led to 
wide acceptance of the model with an accreting black hole 
in the nucleus, and it offered an unprecedented opportunity 
to explore directly the pattern of the accretion flow onto the 
central hole in active galactic nuclei (AGN; for recent and 
detailed expositions of the subject, see Peterson 1997; Krolik 
1999). The iron reflection features and remarkable variabil- 
ity patterns have been reported in many other AGN, galactic 
black-hole candidates, and similar objects. 

It is the main aim of the present paper to examine re- 
lationship between the form of motion of the gaseous ma- 
terial in a galactic nucleus and the resulting proflle of the 
spectral features. We concentrate ourselves on the question 
whether current observational evidence can be explained in 
terms of individual clouds with spherical or almost spherical 
(rather than disc-type) orbital motion around the centre, or 
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if strongly flattened distribution of the clouds is preferred re- 
sembling a disc or a ring (different models of the line forma- 
tion in AGN were reviewed by Netzer 1990). This topic has 
far-reaching consequences for the unification scheme of AGN 
(Antonucci 1993). Although the assumptions of the disc ge- 
ometry and of strictly planar bulk motion of the gaseous 
material are relaxed in the present paper, the presence of a 
compact supermassive accreting nucleus remains crucial for 
explaining the observed spectral features. 

Broad iron features were expected in X-ray spectra on 
the basis of the model in which the iron line is formed on 
the surface of a geometrically thin, optically thick and rel- 
atively cold medium after irradiation by a primary source 
(Fabian et al. 1989). Subsequent detailed fits to observa- 
tional data confirmed that the line profile is in agreement 
with the predictions of the accretion disc model in which 
the intrinsically narrow line, emitted at energy 6.4 keV, is 
gravitationally shifted and Doppler broadened/boosted due 
to the disc orbital motion. Comparisons of the model with 
the data allow one to constrain the disc inclination angle, 
the range of disc radii contributing to the emission, and 
the radial dependence of the incident X-ray flux, although 
there are various uncertainties if astrophysically more re- 
alistic models of accretion flows are introduced, and if the 
lack of resolution and the noise in available data are taken 
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into account; cf. Fabian et al. (1989); Tanaka et al. (1995); 
Weaver & Reynolds (1998) for the case of iron line diag- 
nostics, and Rokaki & Boisson (1999) for application to UV 
continuum and H/3 emission line. 

The expected amplitude of the line profiles and the 
characteristic form of continuum from X-ray illuminated ac- 
cretion flows were examined, taking into account combina- 
tion of effects due to high orbital velocities and strong grav- 
ity near the nucleus (George & Fabian 1991; Matt, Perola & 
Stella 1993). It has been argued with various levels of refine- 
ment that parameters of the central black hole (especially its 
angular momentum) can be inferred from disc-line spectra, 
assuming that they are sensitive to the radius of the inner- 
most stable orbit whose imprint is visible in radiation of the 
accretion flow (Laor 1991; Iwasawa et al. 1996; Dabrowski 
et al. 1997; Pariev & Bromley 1998). 

Reynolds & Begelman (1997) pointed out, however, that 
the model is not that sensitive to the value of the spin of the 
hole if the contribution to the line is allowed also from mat- 
ter inspiralling below the marginally stable orbit, rms- Such 
a possibility appears as perfectly consistent with high effi- 
ciency of X-ray radiation if the adopted geometrical depth of 
the stream is small, in agreement with actual computations 
of the flow properties below rms (Muchotrzeb & Paczynski 
1982). 

In the case of accretion-disc geometry, and with differ- 
ent assumptions about the X-ray illumination and repro- 
cessing, the predicted spectra (line plus continuum) have 
been calculated by several authors (recently Young, Ross & 
Fabian 1998; Martocchia, Karas & Matt 2000; see further 
references therein). The same scheme, coupling the shape 
and the amplitude of the line with the shape and the ampli- 
tude of the Compton reflection component, was successfully 
applied to galactic black-hole candidates (Zycki, Done & 
Smith 1998; Done & Zycki 1999). 

The line flts are consistent with the line emitted at in- 
trinsic energy 6.4 keV (so it comes from weakly ionized iron) 
while the broad shape of the spectral feature indicates that 
it is formed near the innermost part of the disc, following in- 
tense irradiation. This might be related to the specific slope 
of the hard X-ray emission (Rozaiiska et al. 2000; Nayakshin, 
Kazanas & Kallman 1999). The blue wing of the line is in all 
sources linked with the broad feature around 6.4 keV, requir- 
ing somewhat special interplay between model parameters. 
Namely, strong constraints are imposed on the inclination 
angle of the disc [see Guainazzi et al. (1999) for the case of 
MCG-6-30-15; Wang, Zhou & Wang (1999) for NGC4151; 
Nandra et al. (1999) for NGC3516] although contribution 
to the total light from various components helps to ease this 
constraint to some extent. 

Sulentic et al. (1998b) raised the problem of disagree- 
ment between the inclination angle derived by two indepen- 
dent approaches: from the iron-line model, and from H/3 and 
Ha measurements. In a large sample of objects, the position 
of the line centroid at 6.4 keV is not consistent with the 
random orientation of the disc with respect to the observer 
(Sulentic, Marziani & Calvani 1998a) although this prob- 
lem is weakened by the fact that the Balmer lines do not 
necessarily have to come from the outer parts of the disc, 
as assumed in the paper. In particular, the mean inclination 
derived for a sample of Seyfert 2 galaxies does not differ from 
the mean inclination angle of Seyfert 1 galaxies (Turner et 



al. 1998), which disagrees with the widely preferred unifica- 
tion scheme (however, see Weaver & Reynolds 1998). Also, 
the outer disc radius comes out lowish, of the order of ten 
Schwarzschild radii (rg=2.95x 10^ M/Mq cm in terms of the 
central black-hole mass M). The problem of this apparent 
disagreement in inclination angles can be weakened by the 
fact that the Balmer lines do not necessarily have to come 
from the outer parts of the disc, as assumed in the paper, 
and there is also a contribution from the narrow unresolved 
component (expected to arise in the dusty /molecular torus; 
Krolik, Madau & Zycki 1994) which has not been properly 
acounted. 

Certain doubts concerning the disc model for the iron- 
line production revived the interest in the alternative expla- 
nation of the line profile (Czerny, Zbyszewska & Raine 1991). 
Misra & Sutaria (1999) assume that the line is produced by 
Compton scattering of the line photons in warm, Thomson 
thick material. They show that such a model fits the data 
equally well as the disc model, but the approach may look 
rather ad hoc because the origin and location of the Comp- 
tonizing medium have not been addressed in their paper, 
neither the source of the X/UV incident continuum that is 
a necessary ingredient of this model. A better motivated ap- 
proach was adopted by Abrassart (2000a; see further refer- 
ences cited therein) who explored spectral properties within 
the frame of the clouds model of accretion onto a black hole 
where both the line formation and Comptonization occur in 
the same medium. However, his results did not incorporate 
the global effects (gravitational redshift, clouds motion, etc). 
Hereafter we will argue that these effects cannot be ignored. 

In the present paper we discuss the possibility of ex- 
plaining the Fe Ka line within the frame of the model in 
which the innermost part of the disc is disrupted due to 
disc instabilities. Part of the disc material forms optically 
thick cold clouds, while another fraction heats up to high 
temperatures acting as a source of X-rays. The clouds are 
not confined to the disc equatorial plane, and they form a 
layer covering a significant portion of the sky from the point 
of view of the central X-ray source (CoUin-Souffrin et al. 
1996). In this model the line profiles are determined by two 
concurrent effects: they arise partially from the Comptoniza- 
tion within the material significantly ionized at the surface 
of the clouds (the intrinsic line profile), and partially from 
kinematics of the clouds distribution and from strength of 
the gravitational field where the clouds persist (smearing of 
the profile). In the next section we formulate a simplifled 
model which captures the essence of the clouds scenario and 
can be further developed to a more realistic form. We show 
the predicted line profiles as a function of model parame- 
ters (specified by the clouds distribution) and the observer 
view angle. The intrinsic shape of the spectral feature has 
been taken either as a narrow delta-type line, or a numer- 
ically computed broad feature (corresponding to high ion- 
ization parameter ^). Then we briefly discuss how various 
complications (obscuration of the clouds, non-Keplerian or- 
bital motion) are reflected in resulting proflles (Sec. ^). We 
conclude the paper by comparison with the ASCA data for 
MCG-6-30-15, and we discuss overall advantages and the 
problems of this picture. 
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Figure 1. A schematic view siiowing the main ingredients of the 
modeL Most of the clouds are distributed at some characteristic 
radius, their origin being presumably in the accretion disc dis- 
rupted due to instabilities. The clouds survive as individual enti- 
ties orbiting around the massive and compact core. They are ir- 
radiated from the centre, and the primary X-rays are reprocessed 
on the clouds surfaces. From the observer's view, dark parts of 
the clouds are in front of the nucleus while radiating surfaces cor- 
respond to the clouds in the back side. The intrinsic spectrum is 
determined mainly by the Comptonization, depending sensitively 
on ionization parameter. Rather clean line of sight to the central 
source emerges among the clouds. The observed spectral features 
are then influenced by the clouds distribution, their velocities, 
and self-obscuration effects. Also the aspect angle of the observer 
must be considered if the distribution is flattened rather than 
spherically symmetric. 



2 THE MODEL 

We assume that the iron Ka line is produced by irradiated 
surfaces of the clouds. The clouds could be formed from the 
inner accretion disc which is eroded at a distance of a few or 
a few tens of Schwarzschild radii. The basic principles of the 
model were outlined by CoUin-Souffrin et al. (1996), and the 
broad-band spectra following from this model were studied 
by Czerny & Dumont (1998). The expected shape of Ka line 
was calculated with careful consideration of X-ray reprocess- 
ing (Abrassart 2000b) but without taking into account any 
kinematical and gravitational effects which must influence 
the observed spectral features substantially if the clouds are 
to be located so close to the central black hole. Resulting 
width and centroid energy of the line were thus entirely due 
to the Comptonization. In those papers, the clouds distri- 
bution was assumed to be spherical and positions of all the 
clouds were fixed at a given distance from the black hole. 

Here, we relax the assumption of stationary clouds 
forming a spherical layer with strictly uniform distribution. 
Instead, the clouds follow planar orbits departing from the 
equatorial plane of the disc (Figure ^ . One of the appealing 
advantages of this scheme is the naturally modeled variabil- 
ity of X-ray sources due to obscuration events (Abrassart & 
Czerny 2000). In the present paper it is assumed that ex- 
cursions of the clouds from equatorial plane are limited by 
the maximum angle of deflection, say Om. This angle, and 
radius of the clouds belt, r, stand as two free parameters in 
the simplest version of the model. The assumed distribution 
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Figure 2. A representation of the clouds distribution (positions 
and velocities of the clouds). In this sketch we consider all the 
clouds being located on the surface of a sphere. Their trajecto- 
ries are characterized by constant value of the radius r, and by 
maximum excursion from the equatorial plane. The clouds are 
distributed randomly within the allowed range of latitude which 
also restricts directions of the clouds velocities (indicated by short 
dashes). Observer's inclination angle 6o is also shown. 

accommodates both the case developed by CoUin-Souffrin 
et al. (1996), and the case of a narrow ring (Gerbal & Pelat 
1981) in the limiting values of Q-^. In principle, our model 
can accomodate a system of clouds on both prograde and 
retrograde orbits with respect to the disc, but an intermedi- 
ate value of the clouds departure from the equatorial plane 
appears as most favourable, corresponding to quasi-spherical 
or somewhat flattened system of the clouds following more 
or less corotating trajectories. This case allows for quite high 
covering factor (necessary in the model which relies on multi- 
ple scattering of radiation among the clouds) while reducing 
the rate of collisions (which otherwise tend to destroy the 
the clouds). 

The present model does not require all the individual 
clouds to preserve their identity for many orbital periods 
or even indefinitely. Although it may appear easier in this 
scheme to think of the clouds as separate entities, the pre- 
dicted spectrum will not change if the clouds have a finite 
and very limited lifetime, assuming they are continously re- 
placed by newly created ones. Each of the clouds can be 
removed from the system after, say, one revolution, either 
due to coUisional destruction with other clouds or by heat- 
ing up to the state when it does not contribute to the line 
spectrum and possibly evaporates. In the other words, the 
spectral features would be expected more pronounced dur- 
ing the initial phase of the cloud's history before reaching 
the equilibrium. 

Radiation of the clouds surface is characterized by in- 
trinsic emissivity which is assumed constant and isotropi- 
cal in the local comoving (rest) frame of each cloud. Since 
the distribution is not spherically symmetric in general, the 
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observed spectrum depends on inclination angle of the ob- 
server. We neglect the contribution to the line coming pos- 
sibly from the outer, fairly smooth part of the disc. The 
intrinsic line profiles of the clouds come out somewhat dif- 
ferent from those which are usually adopted in the disc-line 
model (cf. Sec. and the random motion of the clouds 
must be also taken into account in calculations of predicted 
spectra. 

2.1 Kinematical and gravitational effects 

The basic observational properties relevant for this paper, 
i.e. the observed profiles of spectral features (their centroid 
energies and widths in particular) , can be roughly predicted 
by a simplified model in which the clouds distribution is 
specified by a small number of phenomenological parame- 
ters. Notice, however, that one still needs to solve for the 
radiation reprocessing in order to obtain the intrinsic line 
profiles. 

In the simplest version of the model, the observer is lo- 
cated at inclination 60 far from the source (0 < 60 < 7r/2), 
while all the clouds are distributed at the same radial dis- 
tance and they move with the corresponding Keplerian or- 
bital velocity, «k(r). Gravitational field of the central body 
is spherically symmetric and described by the Newton law. 
We introduce Schwarzschild radius = 2GM/(? as a 
convenient length-scale and we use geometrized units with 
G = c = 1 hereafter. 

The clouds are distributed randomly within a band of 
spherical latitudes: 7r/2 — < 8 < 7r/2 -I- (Figure |^. 
Spherical coordinates r, 6, (j) ''■re employed with the disc 
plane a.t 6 = -n /2. Parameter 6'ni refiects the dynamics of the 
process of clouds formation and of their interaction with 
surrounding environment. If the disc disruption proceeds 
rapidly, and if the newly formed clouds are efficiently cap- 
tured by some outflowing plasma and/or accelerated by the 
radiation pressure, then the clouds distribution can be ap- 
proximately characterized by r = const with relatively large 
6m. On the other hand, if the clouds acceleration is only 
moderate in the direction perpendicular to the disc plane, 
and if their subsequent evaporation is fast in comparison 
with the timescale of latitudinal motion, then the distribu- 
tion is again well characterized by constant radius but rather 
small 6'm. 

Velocity of the clouds points along r = const surface 
and its direction deflnes the angle a: 

« = Uk(cosa 60 — sina ee). (1) 

The range of possible values of a is determined by the an- 
gular width 9jn of the clouds distribution. Straightforward 
trigonometry gives relation for the maximum a of the clouds 
velocities, 

sin a = ± (cos 9 cot 9 cot 9m cos 6m — sin 6 sin 6m ) 
X (^1 — cos^ 6 sin^^ 9m) 

=p cos^ 6 smT^ 6'm (1 — cot^ 9 cot^ 6'm) ^''^ . (2) 

This relation is shown in Figure ^ where 9m stands as pa- 
rameter of the curves a{6\ 6m)- For example: one can read in 
the graph that latitudinal motion of the clouds has turning 
points at a(7r/2 — ^m; ^m) = where trajectories are paral- 
lel to the disc plane (i.e. v \\ e^). The upper/lower signs in 
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Figure 3. The range of allowed values of angle a constraining the 
velocity direction of the clouds as a function of Q [cf. eq. (^]. The 
curves are labeled with the corresponding latitudinal boundary 
dm of the clouds belt, which adopts values from dm = t/S (strictly 
spherical distribution) down to dm = (the equatorial ring). 

eq. (2) correspond to the descending/ascending parts of the 
trajectory, respectively. We recall that for 0m ^ the clouds 
remain in the equatorial plane (ring-like distribution) , while 
for 9m = 7r/2 the clouds are spread uniformly over the whole 
sphere (with random distribution). 

The observed radiation of a cloud is determined by 
special-relativity effects which influence photon energy 
E and intensity 1(E) in usual manner (Cunningham & 
Bardeen 1973): I{E) = h{E / g) where h is the emitted 
intensity, and g is the redshift factor deflning the change of 
energy. We adopt standard notation, similar to Ohna et al. 
(1995) who studied a related problem with application to 
fast winds. In our case, however, it is the inner part of the 
clouds surfaces which is hot enough to emit the X-rays. The 
remote part of the clouds (with respect to the central source) 
is not subject to primary irradiation, and it emits almost no 
X-rays at all. As a consequence of self-obscuration. X-rays 
are visible only from the clouds with appropriate orientation, 
and their radiating area is further reduced on the rim of the 
image by projection effect. In addition, mutual eclipses of 
the clouds must be considered. The resulting obscuration is 
enhanced in the case of clouds which substantially overlap 
each other; we will characterize this effect by another free 
parameter [ui; cf. eq. (^ below]. 

2.1.1 The case of negligible obscuration 

First we ignore entirely the obscuration of the clouds by 
other clouds. We assume that the clouds radiate isotropi- 
cally in their local comoving frame, and they are of spheri- 
cal shape with an inflnitesimally small radius (i.e. much less 
than Tg). In this case all the clouds from the whole hemi- 
sphere contribute to observed X-rays. 

The redshift factor g depends on the angle between ve- 
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locity of the cloud and the unit vector s towards the ob- 
server, 

g^L-\{l^v.s). (3) 

Here, 7 = — is the Lorentz factor, and v.s — 

Wk cos 9o sin a sin 6* — iik sin 60 (cos a sin + sin a cos cos 0) . 
Integration of hght over the clouds distribution yields the 
total observed radiation flux. Effects of general relativ- 
ity are acc ounted by the gravitational redshift term L = 
•^1 — r^/r, but gravitational lensing has been ignored here 
for simplification (this might have some effect on the light 
from clouds at the upper conjunction). 

Since we deal with the clouds that are irradiated by the 
central source, one expects the local emissivity to be depend- 
ing on the angle between the direction towards center and 
the normal to the cloud's surface. In such a case, by a simple 
geometrical argument, the observed line proflle of the whole 
axisymmetric collection of the clouds is the same as for the 
case of clouds radiating truly isotropically. This conclusion 
does not hold, however, if partial and non-axisymmetric self- 
obscuration of the clouds is involved. 

2.1.2 Partial obscuration of the clouds 

The presence of a large number of clouds (typically A'' ~ 10^) 
enhances the chance of mutual obscuration among them. 
This effect cannot be ignored in this model because the ob- 
server receives X-rays only from those parts of the clouds 
which are on the inner side of the r = const surface, and 
not too close to the edge of the image (as projected onto 
the observer image plane). It can be included in the model 
by introducing another parameter, say ui, characterizing the 
maximum angle between the line of sight and the edge of the 
clouds distribution: s.e^ = cos(7r — ui). In other words, out- 
going rays must pass from interior hot surfaces of the clouds 
through empty holes in between them before they can fi- 
nally escape towards the observer. This is impossible near 
the projected rim of the sphere. Therefore in calculations 
of the observed radiation, only those clouds are considered 
which, in addition to the conditions described in the previ- 
ous paragraph, satisfy also relation (Figure ^) 

sin 9o sin 9 cos -I- cos 9o cos 9 > — cos oj. (4) 

Let us remark that the angle uj can be related to the 
covering factor /c (defined as the fractional area of the sky 
subtended by the clouds as viewed from the centre of the 
source). The covering factor /c can be expressed in terms 
of the clouds typical diameter d, typical distance I between 
the clouds, and their number A'^: /c = N{d/r)'^, with the 
obscuration constraint d~l cos u} for spherical clouds filling 
the whole r = const surface (/c obviously decreases with 
9in decreasing). The exact description of obscuration would 
require a specification of the clouds shapes, and it can be 
carried out only numerically by generating random clouds 
distributions with the constraints imposed by the visibility 
conditions, checking the line of sight in each case. Further- 
more, d and / are linked to each other via A'^, but parameter lj 
has direct geometrical meaning in our model, and it appears 
thus convenient for the present discussion. In principle, the 
same value of the phenomenological parameter ui can corre- 
spond to different physical models of the clouds origin, their 
geometrical configuration and other details. 




Figure 4. A view of the clouds distribution in the case of non- 
negligible obscuration. Positions and velocities of those clouds are 
shown which fulfil conditions of visibility (^. In comparison with 
Fig. ^, only a fraction of the clouds contribute to the observed 
spectrum: those moving more or less perpendicularly to the ob- 
server's line of sight. Direction towards the observer is shown with 
inclination angle 60. Different regions are indicated on the sphere 
according to visibility of the clouds. The clouds do not reach high 
latitudes, and they enter in dead angle near the rim of the pro- 
jected image. This figure corresponds to parameters dm = 30°, 
60 = 60°, u! = 60° in our notation. 

2.2 The intrinsic line profile 

The intrinsic shape of the iron line has been computed using 
a Monte Carlo code noar (Abrassart 2000b), coupled with 
the code titan (Dumont, Abrassart & Collin 2000) com- 
puting the ionization state of the gas for the assumed value 
of ionization parameter ^=F^/n {F-^ denotes the incident 
X-ray fiux, n is average number density of a cloud). The 
code TITAN serves to compute the radiative transfer in the 
Compton thick medium in a broad-band frequency range, 
and it includes the computations of opacity and tempera- 
ture structure. It also gives local emissivity of the gas in 
various components of the Ka line. The NOAR code uses the 
resulting opacity, it computes, in more detail, the transfer 
of hard X-ray photons, and provides better description of 
the details of hard X-ray spectra including Comptonization 
of the line within the hot surface layers of the clouds. It 
also describes the energy deposit by hard X-ray Compton 
heating for titan. The two codes are thus coupled and they 
both together offer the most accurate description of the iron 
line emission from the irradiated medium, assuming it is 
optically thick for electron scattering. 
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Figure 5. The intrinsic emissivity profile (arbitrary units) around 
the 6—7 keV spectral feature, as predicted by the Comptonization 
model for hot parts of the clouds surfaces. The three curves corre- 
spond to different values of the ionization parameter ^ (Abrassart 
2000b). 
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Figure 6. Relation between radius (logarithmic scale) and cen- 
troid energy of the line features dominating the intrinsic spec- 
trurn, after multiplication by that part of the redshift factor from 
eq. (H) which depends on r only and causes the overall redshift 
of the line. These curves correspond to the same values of § as in 
Fig. I 
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Figure 7. The expected line profile produced by the clouds 
sphere with r = 8M. The intrinsic emissivity has been assumed 
as a delta line at 6.4 keV. The four panels correspond to different 
angular widths of the clouds belt: Bra = 5° (almost equatorial 
ring), 30°, 60°, and 90° (strictly spherical distribution). In each 
panel, four curves correspond to different values of the observer 
inclination: 9o = (pole-on view, therefore the narrowest pro- 
file), 30°, 60°, 90° (equator-on view, the broadest profile). No 
obscuration is involved here. See the text for details. 



line comes out almost unresolved, and any smaller degree of 
ionization would result in even narrower line. The line for 
^ — 10"' shows the peak shifted to 6.7 keV, and a broad red 
wing. The line for ^ = 10* contains a strong contribution of 
6.9 keV component arising from hydrogen-like atoms; its red 
wing is extremely broad for such a high value of ionization. 

Since the numerical computations of intrinsic profiles 
are exceedingly time consuming, we detach the line from its 
underlying reflected continuum and discuss the line spec- 
trum separately. Such a simplification is possible because 
the line is not very sensitive to the slope of the incident X- 
ray radiation, while the reflected continuum depends on the 
slope considerably and would require computations of a set 
of different incident continuum slopes. 

Figure ^ shows how the intrinsic centroid energy of the 
line is moved towards lower energy by factor L{r) 'y~^{v\i{r)) 
in eq. (^). This term is function of radius only. Therefore 
the graph shows roughly the mean energy of the expected 
observed line peaks from Fig. ^ depending on the radius of 
the clouds sphere. 

In the next section we calculate the expected profiles 
for different parameters of the model. In particular, we fur- 
ther discuss the effect of broadening of the line, which arises 
from (;6-dependent term v.s in eq. (^ and results in large 
equivalent widths and double-horn profiles in some cases. 



The computations were performed for a plane parallel 
slab corresponding to a sphere with large radius, so that 
multiple reflections between the clouds were neglected (they 
can effectively change only the value of the ionization pa- 
rameter). The line profiles used in our computations of the 
cloud model are shown in Figure ^ The pattern obtained for 
^ = 10^ is centered at 6.4 keV, and it is so narrow that the 



3 RESULTS 

3.1 A narrow intrinsic line 

The computations performed with the codes titan and 
NOAR gave the Fe Ka line strong but practically unresolved 
for the case of ionization parameter f = 100 under the 



© 0000 RAS, MNRAS 000, 1-12 



A cloud model of active galactic nuclei 7 



r/M = 20 






0.4 7 

Energy 



7 10 
Energy 



Figure 8. The same as in Fig. ^ but for r = 20M. Sucii intrinsi- 
cally narrow profiles are expected in the case of low ionization of 
the clouds. They can also serve as templates for convolution with 
broader and more complex intrinsic profiles from a highly ionized 
medium. 
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Figure 10. The predicted spectral profiles corresponding to 
delta-type line as in Fig. but now the obscuration parameter 
has been introduced with w = 30° . The dependence of the profiles 
on observer's inclination is substantially reduced in comparison 
with the previous case of negligible obscuration. Here, r = 8M. 
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Figure 9. The same as in Fig. I but for r = lOOM. 
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Figure 11. The same as in Fig. but for r = 20M. 



adopted energy grid: AE ^ 0.07 keV, and AE/E ^ 0.01. 
A faint red wing is barely visible and it contains less than 
a few per cent of the line flux. Therefore, in order to study 
the case of a neutral iron line we could neglect the intrinsic 
width of the line and we performed computations assuming 
that all photons were emitted with energy 6.4 keV (delta- 
type line with no continuum background) . The observed line 
profile is thus determined exclusively by radius of the clouds 
distribution r, the level of departure of the clouds from the 
equatorial plane Om, observer's angle 9o, and the obscuration 
parameter u which is further related to the covering factor. 

We show the predicted profiles for three different radii 
r of the clouds distributions (Figures |^-^. If the clouds 
occupy a very narrow belt on the sphere (^m ^ 5°), the 
resulting profiles are reduced to those obtained for a single 



ring within the frame of the disc model of the line formation 
(Gerbal & Pelat 1981; Laor 1991). The corresponding values 
of r and 9^ are given on top of each figure. Two horns of the 
line are clearly seen with unequal heights due to the Doppler 
effect. The dependence of the profiles on 9o gets gradually 
diminished with increasing sphericity of the clouds distri- 
bution, as can be checked by comparing four panels with 
different 9ni- The profiles are normalized to the maximum 
flux. 

Figs. correspond to zero obscuration of the clouds 
{lu — 90°); the clouds do not shield each other and they 
all are seen by the observer. We remark that the resulting 
profiles are not completely smooth in these graphs because 
they are produced by individual sources of light. Naturally, 
the curves would appear very wiggly if the number of clouds 
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and the corresponding covering factor were small, but our 
simulations indicate quite large covering factors. 

The resulting line profiles come out substantially differ- 
ent from the previous case of zero obscuration if uj <^ n /2: 
the profiles get narrower, and they loose, partly or even 
completely, their double-horn shapes, so characteristic for 
ring- type sources (Figures p^-pj]). Considered jointly with 
the underlying continuum, there is less power in the line, so 
that equivalent width is diminished. The dependence on ob- 
server's inclination is further reduced by obscuration. One 
can easily deduce that the profiles do not depend on Oq if uo 
is small enough: uj < 9q -\- 9^ — 7r/2. Also, notice that the 
constraints on the clouds visibility are not satisfied in case 
of very small ^o, ^ and 9m, so that the curves are shown 
only for appropriate combinations of these parameters. 

What remains for the next section is to complete the 
above discussion by computing the expected profiles of more 
complicated intrinsic emissivities. 

3.2 A broad intrinsic line 

Now we examine the spectral features with intrinsically 
broad profiles and high ionization parameter ^ 10''; 
Abrassart 2000b). In each case the continuum had been fit- 
ted by a power-law outside the feature (4-9 keV) , and then 
the result was subtracted from the predicted spectrum. In 
this way the expected spectral features have been obtained 
(Figures p2|-p^). We recall that the full (line and contin- 
uum) intrinsic spectra are the result of radiation transfer 
around 6.4 keV with multiple scattering among the clouds; 
they differ substantially from those corresponding to simple 
delta- line profiles. Here we present narrow band spectra; the 
actual energy range is shown on abscissae. The two figures 
differ by obscuration parameter u): obscuration was ignored 
in the former, while it was assumed to be very large in the 
latter case. 

Let us note that the values of /c inferred from spectra 
fitting procedures (and corresponding parameter ua of the 
model) represent limiting boundary estimates of these quan- 
tities. One of the reasons of uncertainty is the clouds velocity 
which until now has been assumed fixed and equal to v-^{r) 
while, in case of strong radiation pressure, the clouds mo- 
tion may be non-Keplerian. The layer of clouds can be par- 
tially supported by intercepted radiation flux coming mainly 
from the center, so that steady-state motion of the clouds 
becomes sub-Keplerian. Therefore, given the observed line 
width, the actual covering factor will be less than that re- 
quired on the basis of clouds free orbital motion. In the case 
of quasi-spherical distribution, the effect of radiation pres- 
sure is maximal for Eddington-type equilibrium condition, 
GMnip — A, where nip is the proton mass, and 



is the outward-directed radiation support due to total lumi- 
nosity L corrected for partial obscuration. For large covering 
factor the total luminosity is small (/c —> 1, L ^ 0) and the 
ratio on the right-hand side of eq. (^ remains finite. Cen- 
troid energy of the line from such a static system of the 
clouds is determined by the intrinsic spectrum and gravita- 
tional redshift only, and the observed line width is reduced 
to the intrinsic one. 



r/M = 20 




Energy Energy 

Figure 12. The expected iron-line profiles (power-law continuum 
subtracted) produced by the clouds at r = 20M for four values of 
the inclination angle of observation 9o (0°, 30°, 60°, 90°). Again, 
the panels are labeled with the adopted values of the width of 
the belt, 6^- The angle w = 90° corresponds to zero obscuration 
factor. The intrinsic spectrum is also plotted with the maximum 
at 6.7keV, corresponding to the assumed ionization parameter, 
? = 103. 



r/M . 20 




3.3 The case of MCG-6-30-15 

The Seyfert 1 galaxy MCG-6-30-15 is a notable example 
in which the extensive observational material has been col- 
lected. It shows the primary spectral slope of F ~ 2.00±0.05, 
as determined on the basis of joint ASCA/KJCTE observa- 
tions lasting « 400 ksec with the cutoff energy at ~ 100 keV 
(Lee et al. 1999). The disc-line scheme has been rather suc- 
cessful also for this object (Tanaka et al. 1995) and the con- 
clusions drawn from the ASCA data are in agreement with 
the analysis based on BeppoSAX observations (Guainazzi et 
al. 1999). As a consequence, strong limits can be imposed on 
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Figure 14. The continuum and the Fe Ko line profile from the 
model cloudfe fitted to the ASCA data for MCG-6-30-15. Re- 
sulting parameter values are r = 29.1, T = 1.856. The remaining 
parameters were kept fixed: ^ = 10*, dm = 90°, = 0.1. 




channel energy (keV) 

Figure 15. A detail of the line and the best fit profile are shown 
with better resolution between 4—8 keV. 



alternative models (Fabian et al. 1995; Reynolds & Wilms 
2000). However, in view of the complexity of the problem, 
it still seems somewhat premature to condemn other mod- 
els as nonviable. We also remark that this object is probably 
somewhat exceptional case, though it is the best studied one 
at present. 

We tested our model of the line profile predicted by the 



clouds model by fitting it to the ASCA data from July 1994 
(Tanaka, Inoue & Holt 1994; Fabian et al. 1994; Tanaka et 
al. 1995; Ebisawa et al. 1996). In order to avoid problems 
with modelling of the reflected component, the contribution 
from the warm absorber (Otani et al. 1996), and the cold 
absorption along the line of sight, we restrict our discussion 
to the SIS data in 3-10 keV band where the continuum can 
be represented roughly by a single power law. Since the de- 
tailed fitting is beyond the focus of the present work, we do 
not consider various luminosity states separately (instead, 
we follow the initial approach of Tanaka et al. 1995) . 

In order to have a reference for the results from our 
model we started by fitting the ASCA spectrum with a 
power law and a disc-line component (Fabian et al. 1989) 
from XSPEC (Arnaud 1996); see Table for the results of 
statistics. We assumed the line energy at 6.4 keV, con- 
sistently with the previous results. The X-ray incident fiux 
was allowed to vary with the distance as ocr"'' with inner 
and outer disc radii fixed at 6 GM/c? and lQr'GM/(?' respec- 
tively. Then we compared those results with the fits to our 
model in two different parameter ranges. 



3.3.1 The spherical configuration 

First we consider the case of spherically symmetric arrange- 
ment of the clouds. The clouds cover the entire sphere 
{9m = 90°), so that the inclination angle is of no impor- 
tance. Fitted parameters are: the photon index of the power 
law r, radius of the clouds distribution r, and the visibility 
factor vj defined as 



V/2' 



(6) 



Different versions of the model cloudfe were examined 
using XSPEC supplemented by three local model subroutines. 
They differ by corresponding values of ionization parameter: 

(i) an intrinsically narrow line in the rest frame of the 
clouds; this case refers to low ionization of the clouds sur- 
faces where X-rays are reprocessed, ^ ^ 10^; 

(ii) a broad line with ^ — 10"^ of the local emission; 

(iii) a broad line with ^ — 10* of the local emission. 

The first of our models, with an intrinsically narrow 
line, does not provide adequate representation of the data. 
Strong residuals remain (as in the case of the Gaussian 
model) , and the red wing cannot be reproduced successfully. 
However, for intrinsically broad line (with 10'^ ^ ^ ^ 10*) 
the red wing is formed locally as the result of Comptoniza- 
tion, and it has roughly the required shape. This model re- 
sults in comparable to the diskline profile (Tab. 
The larger value of the ionization parameter gives marginally 
better fit than the lower one. Results of the fitting procedure 
with SISO and SISl detectors are shown in Figure JlJ. The 
line is reproduced with better resolution in Figure (here 
we show only SISO data for clarity). 

Our calculation imposes limits both on radius of the 
clouds distribution and on their covering factor. The for- 
mer is constrained by requirement of appropriate location 
of the peak of the line; the gravitational shift has to bal- 
ance somewhat increased energy of the line emission in this 
model (Fig. ^). The latter dependence reflects the fact that 
small radius of the clouds distribution goes in hand with 
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Table 1. Spectral fitting of the model parameters to the data for MCG-6-30-15. Standard statistics was examined for different XSPEC 
models. The first two lines in the table refer to standard models for the sake of comparison (diskline and powerlaw), while the rest 
corresponds to the user supplemented model (cloudfe) as described in the text. The resulting values of three fitted parameters are given. 



Model ^ Parameter #1 Parameter #2 Parameter #3 x^/d.o.f. 



diskline 




r = 


1.8iU_o 


o _ _2 4t;+0.34 


eo = 31.5j;3;2deg 


723.6/720 


powerlaw 




r = 


1 79O+0-020 
^- '^'^-0.021 






852.5/722 


cloudfe, spherical 




r = 


1 700+0.021 


r = 38tJ*i GM/c2 


ro = 0.52^0 39 rad 


804.1/720 


cloudfe, spherical 


10^ 


r = 


1 702+0-024 
1. 022 


r = A8Atl°^^ GM/c^ 


ro < 0.19 rad 


725.3/720 


cloudfe, ring-type 


10^ 


r = 


1 796+0025 


r = AQ.Otf.f GMjc^ 


< 21.7 deg 


724.9/720 


cloudfe, spherical 


10* 


r = 


1 Sc;fi+0 029 

1.00D_g Q29 


r = 29.lti-lGM/c'^ 


ro < 0.13 rad 


720.0/720 


cloudfe, ring-type 


10* 


r = 


1 oc:7+0.026 
I.c50/_Q 024 


r = 29.8+4 GM/c2 


6»ni < 19.4 deg 


721.6/720 



dramatically strong smearing of the line, unless the effect of 
self-obscuration prevents us from having a direct view of the 
clouds with the highest velocities along the line of sight. We 
recall that the above discussion assumes Keplerian motion 
of the clouds; numerical values given in the column -oo of 
Tab. ^ must be increased (less obscuration) if sub-Keplerian 
motion of the clouds is taken into account, as discussed in 
previous section. 



3.3.2 The ring-type configuration 

In order to see whether we can detract the strong require- 
ments on the clouds obscuration at the expense of intro- 
ducing partially flattened, ring-type geometry, we fitted the 
model assuming zu = 0.5 and varying 9ni. To illustrate the 
procedure, let us fix inclination angle at 6o = 30°, and 
consider two values of the ionization parameter, ^ — 10^ 
and C = 10*. The results are referred as ring-type models 
cloudfe in Tab. ^. The quality of the fits is comparable to 
those obtained for spherically symmetric case, however, the 
configuration comes out very flat, almost a ring indeed. We 
could not fit the data assuming vj ^ 1 (the case of negligible 
obscuration) because such a flat distribution of the clouds 
clearly requires a substantial range of radii contributing to 
the spectrum (rather than a narrow ring). This conclusion 
is valid independently from the intrinsic broadening of the 
line due to internal Comptonization. 

On the basis of several other attempts to flnd acceptable 
fits with different obscuration, we reached the conclusion 
that the currently available data is still not good enough to 
distinguish definitively in this model between a flat geometry 
and a spherically symmetric geometry of the accretion flow. 



4.1 The role of high ionization 

Since the red wing of the line is mostly due to Comptoniza- 
tion within the clouds surface layers (instead of kinematics 
of the accretion disc), the observed variability of the line 
profile must be also linked with the changes of the physi- 
cal state of the clouds surfaces (Ebisawa et al. 1996; Lee et 
al. 1999). Self-consistent calculations of the resulting spec- 
tra present a difficult task because the thermal stability of 
the gas irradiated by hard X-ray photons depends on many 
factors at temperatures 10^-10® K (Krolik 1999; Rozanska 
1999, and references therein). Even a minor change in the 
incident radiation may cause a dramatic response of the tem- 
perature and, consequently, the density of the irradiated gas 
through its expansion/contraction under constant pressure. 
Such fiuctuations mean a significant change of the ioniza- 
tion parameter and of the intrinsic shape of the emitted 
line, without major alteration of the total luminosity of the 
source or the total fiux in the iron line. 

The demand of high ionization of the medium respon- 
sible for the line emission makes a difference from almost 
neutral environment of illuminated disc-type accretion flows, 
and there is hope to flnd its imprint also in the continuum. 
However, the limited spectral range of ASCA does not allow 
us to exploit such a possibility and to differentiate between 
the disc and the clouds models. 

Also, we cannot presently distinguish the spherical con- 
figuration with strong self-obscuration of the clouds from 
the case of much flatter geometry and less obscuration. A 
better spectral resolution is needed in order to make the 
X-ray spectroscopy really a powerful tool and to determine 
the form of the accretion flow. Further work on the models 
is also needed to demonstrate true differences between the 
models. 



4 DISCUSSION 



4.2 The role of hydrostatic equilibrium 



The effect of obscuration between the clouds plays an im- 
portant role in the model: it restricts line-of-sight velocities 
of the clouds and suppresses the dependence of observed 
line profiles on inclination in comparison with the standard 
formulation of the disc-line model. Another interesting at- 
tribute of the clouds model is the requirement of high ion- 
ization of the material which is responsible for the iron-line 
emission: ^ <; 10''. Such a high ionization has two conse- 
quences which we can briefly outline. 



The calculations of local emissivity (Sect. 2.2) have been 
carried out under constant density approximation. If the 
clouds have enough time to achieve the hydrostatic equilib- 
rium, the radiative transfer should be performed assuming 
constant pressure instead of constant density within the re- 
processing medium. Computations under the constant pres- 
sure assumption or, more generally, under the assumption of 
hydrostatic equilibrium show a specific temperature profile 
due to the radiative thermal instability at the intermedi- 
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ate temperatures (see e.g. Krolik, McKee & Tarter 1981; 
Rozariska & Czerny 1996; Rozariska 1999; Nayakshin et al. 
1999). The irradiated medium is at the inverse-Compton 
temperature in the upper layers while being cold inside, with 
rather sharp transition determined by the effect of heat con- 
duction (Rozahska 1999). Nayakshin et al. (1999) argue that 
the intermediate stable branch at the temperature about 
10® K has a negligible influence on the formation of spectral 
features, while preliminary results of Monte Carlo simula- 
tions (Zycki ct al., in preparation) indicate that a possible 
enhancement of the iron abundance (crucial for the devel- 
opment of this branch; Lee et al. 1999) by a factor of ~ 2 
may change this conclusion. 

The problem remains unsettled for the moment, how- 
ever, if the inverse- Compton temperature itself is low enough 
to allow for the formation of the Ka line, then the line forms 
in a broad zone and the results of the computations do not 
depend crucially on the assumed density profile. This is quite 
possibly the case for the source MCG-6-30-15. 

4.3 The role of inverse-Compton temperature 

The spectral shape in the hard X-ray band was well de- 
termined for this object from the SAX data (Guainazzi et 
al. 1999). The best fit model gave the photon index F = 

2.04 and the high-energy cut-off at 130 kcV. The inverse- 
Compton temperature corresponding to such a spectrum is 
equal to 2.02 x 10^ K, however, the soft photons emitted by 
the clouds and observed as the Big Blue Bump component 
reduce this value further. 

The object is highly reddened in the optical/UV band. 
The extinction estimate based on the Balmer line indicates 
the E(B—V) value in the range of 0.61-1.09 (Reynolds et 
al. 1997). Assuming a minimum value for the extinction 
these authors determine L(NIR— UV)= 0.3 L{X). Somewhat 
larger values arc favoured because too low bolomotric lumi- 
nosity of the internal source makes a large far-IR luminosity 
of this object difficult to explain. The upper limit on the 
extinction rises the flux at log(i/) = 15 by a factor of 12 
in comparison with the lower limit. This fact together with 
an unconstraint high-frequency extension of the Big Blue 
Bump introduce significant uncertainty in the shape of the 
broad-band spectrum. Therefore, we can assume a standard 
shape of the spectrum, as in Abrassaxt & Czerny (2000) , i.e. 
L(NIR-UV)= 3L(X). The observed ratio is slightly differ- 
ent from the intrinsic spectrum seen by the clouds (Abras- 
sart & Czerny 2000; see their Eq. (12)). The correcting fac- 
tor to this ratio is typically equal to 1.05-1.11. Assuming 
the value of 1.1 we finally obtain the Compton temperature 
of 4.7 X 10'^ K as a characteristic value for MCG-6-30-15. 

The derived value falls just within the temperature 
range which resulted from our computations. The upper 
layer of the clouds, down to res = 1, has the temperature of 

6.5 X 10® K assuming the parameter ^ = 10"', and 1.8 x 10® K 
for ^ = 10^. Therefore, in the case of the presented computa- 
tions, the presence or the absence of hydrostatic equilibrium 
does not seem to pose a crucial question. This is quite fortu- 
nate because, under our present insufficient understanding of 
the clouds evolution, it is by no means clear which of the two 
approximations is more appropriate to the clouds that we 
consider - whether it is constant density or rather constant 
pressure assumption. Notice that for other Seyfert galaxies 



than MCG-6-30-15 most probably their inverse-Compton 
temperatures will tend to be higher because their spectra 
are harder. 



5 CONCLUSIONS 

We presented a simple cloud model of AGN in which self- 
consistent computations of the intrinsic X-ray emission are 
supplemented by the phenomenological description of the 
clouds distribution. We demonstrated that satisfactory fits 
to the iron Ka line feature can be obtained, and the model 
should be thus considered as viable and developed further, so 
that the phenomenological parameters are derived from the 
physical model of the clouds origin. We envisage the origin of 
the clouds due to erosion of the inner accretion disc. Detailed 
physics of the clouds formation remains beyond the scope of 
this paper but the adopted parameters (r, 0o, Om, w) should 
refiect this scheme in a natural manner. 

We remark that the scheme described in the present 
model is not restricted by recent considerations of Reynolds 

6 Wilms (2000) because our model falls into that category of 
models in which the X-ray source is not viewed through the 
Comptonizing medium, and the clouds experience somewhat 
different continuum than the observer (Abrassart & Czerny 
2000). 

Observed profiles must be influenced by several other 
effects which were neglected in this paper due to our igno- 
rance of the detailed model of the clouds properties. One 
can however anticipate the expected implications of a more 
refined model. 

First, we ignored gravitational lensing for which one 
can expect a similar conclusion like for partial obscuration 
and sub-Keplerian motion of the clouds: the lensing effect 
enhances radiation coming from the clouds near the line of 
sight and on the remote side of the sphere. If included, the 
effect decreases the covering factor derived here by fitting 
the lino profiles. 

Second, we did not consider radial distribution of the 
clouds. This could be introduced in a straightforward man- 
ner but the assumption of almost constant radius provides 
a very good approximation in the current model. Hot parts 
of the clouds surfaces are those which lie innermost, where 
they are subject to strong primary irradiation by X-rays. 
The impact of the radial distribution can be absorbed in 
parameter w (enhanced obscuration). 

We conclude by emphasizing the fact that the present 
model differs from the disc-line model in the mechanism how 
the intrinsic spectrum is produced. There are however also 
common aspects of the two schemes. Indeed, a ring-like dis- 
tribution of the clouds can be treated as a special case (small 
^m) resulting in profiles comparable to the disc-line model 
with appropriate emissivity. Notice, however, that here we 
obtained reasonable fits under the assumption of spherical 
or almost spherical clouds distribution. Certain degree of 
the fiattening (characterized here by Om) is natural to this 
model where the clouds origin is in continuous erosion of the 
inner disc, but, on the other hand, Om cannot be too small 
because substantial Componization requires large covering 
factor. Finally, we remark that successful fits require strong 
gravitational field of the central body which shifts the line 
centroid towards lower energy. 
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